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A state of an electron in a quantum wire or a thin film becomes metastable, when a static electric 
field is applied perpendicular to the wire direction or the film surface. The state decays via tunneling 
through the created potential barrier. An additionally applied magnetic field, perpendicular to the 
electric field, can increase the tunneling decay rate for many orders of magnitude. This happens, 
when the state in the wire or the film has a velocity perpendicular to the magnetic field. According 
to the cyclotron effect, the velocity rotates under the barrier and becomes more aligned with the 
direction of tunneling. This mechanism can be called cyclotron enhancement of tunneling. 

PACS numbers: 03.65.Xp, 03.65.Sq 



I. INTRODUCTION 

Tunneling in a magnetic field is a matter of investi- 
gation for many years. The magnetic field can influence 
tunneling across a potential barrier in two different ways. 

First, the magnetic field can modify an underbarrier 
motion related to a classically forbidden region. Studies 
of this phenomenon are presented in the literature, fn 
Refs. pj, [2| it was pointed out that an underbarrier fall 
of the wave function can be less rapid in an inhomoge- 
neous sample. See also Refs. [!, HI 0, 111 0] ■ I n ^ e f- H a 
transmission coefficient through a quadratic barrier was 
found. A decay of a metastable state was considered in 
Ref. Q . The certain peculiarities of an underbarrier wave 
function were discussed in Refs. [ToL [TTI] . 

Second, the magnetic field can influence a state of an 
electron at a classically allowed region after an exits from 
under the barrier. A typical example is the Wigner res- 
onance when the electron tunnels into a potential well 
with a level aligned to its energy [12]. See experimental 
measurements, for instance, in Refs. [l3l [LH [l5j . An- 
other example relates to a specific density of states in 
the classical region after the tunneling barrier. A state 
of an electron, influenced by the magnetic field, may fit 
better that density of states and this results in increase 
of tunneling rate [161 ] . 

The goal of the paper is to study tunneling decay rate 
of a metastable state in a magnetic field (the electron 
after tunneling goes to infinity). The question to be an- 
swered: Can a magnetic field increase the decay rate? 
ft is clear that the above effect of density of states at 
the region after the barrier can, in principle, increase the 
rate. But this effect, related to a prefactor, cannot be 
very large. According to Ref. [lfj] , there is approximately 
50% enhamcement of the effect. 

It would be much more amazing to increase the main 
exponential part of the decay rate by the magnetic field. 
The references d S, i, II 1, @, 0, H H say that it is impos- 
sible. Indeed, when an electron enters under the barrier 



its velocity deviates, due to the cyclotron effect, from a 
tunneling path with no magnetic field. This leads to a 
reduction of the tunneling probability by the magnetic 
field. The reduction can be also explained in terms of 
increasing of the total barrier. The additional barrier is 
proportional to a squared velocity of the electron in the 
magnetic field [l|, 0] . 

But there is a situation when the electron tunnels from 
a quantum wire or another object extended in the direc- 
tion perpendicular to tunneling. In this case a state prior 
to tunneling can have a finite velocity perpendicular to 
the tunneling direction. According to the cyclotron ef- 
fect, this velocity rotates under the barrier and becomes 
more aligned with the tunneling direction. This leads to 
enhancement of the tunneling rate by the magnetic field 
(cyclotron enhancement). 

Formally, cyclotron enhancement of tunneling results 
from a reduction of the main tunneling exponent which 
reminds one of Wentzel, Kramers, and Brillouin (WKB). 
The exponent can be reduced in a few times. Suppose 
that at zero magnetic field the tunneling rate is propor- 
tional to 10 -24 . The magnetic field can turn it into, say, 

io- 7 . 

We consider in the paper tunneling from a straight 
quantum wire, directed in the y axis, embedded into a 
two-dimensional electron system in the {x, y} plane. The 
potential barrier is created by the electric field £q di- 
rected along the x axis (the direction of tunneling) . The 
magnetic field H is aligned along z. According to elec- 
trodynamics, after tunneling a motion of the electron in 
perpendicular magnetic and electric fields is restricted 
by a finite interval in the x direction [l7|. To get the 
electron passed to the infinite x one should put some 
potential wall(s) along the x direction restricting the y 
motion. Due to multiple reflections from the restricting 
wall in the magnetic field the electron goes to the infinite 
x. We model the walls by the potential proportional to 

The theory presented relates also to tunneling from a 
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flat {y, z} film with quantized electron motion in the x 
direction. The electron tunnels into a three-dimensional 
reservoir. Restricting walls should be placed parallel to 
the {x, z} plane. 

Without the restricting walls a solution can be ob- 
tained analytically on the bases of the modified WKB 
approach as shown in Sec. lIIII An approximation of clas- 
sical complex trajectories is formulated in Sec. IIVI In 
Sees. Ivl and I VII two different methods of numerical calcu- 
lations are applied to the problem with restricting walls. 



II. FORMULATION OF THE PROBLEM 

We consider an electron localized in the {x,y} plane. 
The static magnetic field H is directed along the z axis. 
Suppose a motion of the electron in the {x, y} plane to 
occur in the potential U(x,y). Then the Schrodinger 
equation, with the vector potential A — {0, Hx, 0}, has 
the form [l2| 



h 2 d 2 ip ti 
2m dx 2 2m \dy 



ix — r— ) ip + U(x, y)ip — Exj), 

where lj c = \e\H/mc is the cyclotron frequency. The 
potential 



At zero magnetic field (h = 0) the WKB probability of 
tunneling [12j | 



w(h — 0) ~ exp 



AB 
"~3~ 



(6) 



is small. In the new variables Schrodinger equation (|TJ) 
has the form 



1 d 2 ^ 
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S(x) 



ip = Aip. 



One can exclude the point x = 0, if to impose the bound- 
ary condition 



dx 



dx 



= -2BV(0,y). 



(8) 



x=-0 



With the semiclassical accuracy, the lowest level at the 
5 well is A = — 1. At a finite electric field that level 
is metastable and we should find a decay rate due to 
tunneling in the magnetic field. 



U(x, y) = —H\ — -5(x) — £ Q x + u 
11 m 



y 



(2) 



describes the quantum wire placed in the y direction 
(the first term), the constant electric field £q (the sec- 
ond term), and the restricting walls in the y direction 
are modeled by the third term. At £q = and H = 
the discrete energy level in the S well (— uq) is a ground 
state in the WKB approximation. 



A. Dimensionless units 

Let us introduce the dimensionless electric field e and 
the magnetic field h by the equations 



u ' 



h = 



(3) 



Below we measure x and y in the units of uq /£q and time 
in the units of 



Too = 



y/2muo 
~£~o~' 



(4) 



The energy is E — uq\ where the dimensionless energy 
A is negative in our problem. We also introduce a large 
semiclassical parameter 



B 



uo^J2muo 
h£o 



(5) 



B. A semiclassical approach 

When the potential barrier is hardly transparent (a 
large B) one can use the semiclassical approximation for 
the wave function [131 



ip(x,y) ~ exp[iBa(x,y)] , 



(9) 



where u(x, y) is the classical action satisfying the equa- 
tion of Hamilton-Jacobi at x ^ 



da\ 2 [da , 
Tx) + {-8-y +hX 



y 



= A. 



(10) 



The function a(x, y) is continuous. As follows from 
Eq. (fT0|) . (da/dx) is also continuous. According to that, 
one can consider the action at positive x with the bound- 
ary condition 



da 

dx 



(11) 



x=0 



which follows from Eq. p|). In the approximation of a 
large B the condition © is reduced to 



dx 



-Bi>(0,y), 



(12) 



x=0 



if to consider the problem at positive x only. 
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FIG. 1: The effective potential energy (I15p is plotted for a 
positive k in Eq. (|16|l . An electron tunnels between x = 
and xo- The two dots mark terminal points between which a 
classical motion occurs. The potential barrier with no mag- 
netic field is shown by the straight line. 



III. CYCLOTRON ENHANCEMENT OF 
TUNNELING 




FIG. 2: (Color online) The decay rate (|22j) of the metastable 
state of the quantum wire versus magnetic field h at B = 15. 
Left curve: The state in the wire has zero velocity parallel to 
the wire, A = — 1. Right curve: The same when the velocity 
is finite, A = -0.594. 



First, we consider tunneling from the quantum wire, 
when there are no restriction of the motion in the y di- 
rection. In other words, we drop down the term y A in 
the potential U(x, y) ([2]). In this case the solution of the 
Schrodinger equation (JTJ) 

ip(x, y) = exp (-iBky) ip(x) (13) 

is determined by the effective Schrodinger equation 

with the effective potential energy 

V{x) = {k-hx) 2 - x -—6(x). (15) 
B 

The first term in Eq. (|15[) is proportional to squared ve- 
locity of the electron. Solutions with positive and neg- 
ative wave vector k are possible according to different 
directions of motion in the wire. In the WKB approxi- 
mation for the state, localized in the S well, 

k = ±Vl + A. (16) 

In the physical units, the velocity along the quantum wire 



= TVl + Aj— . 

m 



2u 



(17) 



Below we consider a positive k related to the motion 
from +oo to — oo in y. For this case the potential energy 
V(x) is shown schematically in Fig.[JJand the WKB wave 
function under the barrier, < x < xq, has the form 



ip(x,y) ~ exp (-iyBVl + X^j 



(18) 



exp 



-B 



dxi 



(x/TTa 



hxi 



Xi 



The classical terminal point, determined by the condition 
V(x) = A, is 



where 



, = W h^ + VTTA) -i. 



(19) 



(20) 



The terminal point xq plays a role of the exit point from 
under the barrier and is marked in Fig. [TJ The classically 
allowed region, V(x) < A, is between two terminal points 
displayed in Fig. [1] Our consideration has a meaning, 
when the classically allowed region is of a finite width. 
This occurs under the condition 



h < 



1 



2(1- Vl + A) 



(21) 



when the two terminal points in Fig. [JJ are separated. 

For the y-independent potential considered the clas- 
sically allowed motion occurs between the two terminal 
points in Fig. [TJ It means that after an exit from un- 
der the barrier the electron motion will remain restricted 
in the x direction. A penetration through the barrier, 
< x < x n , in Fig. [TJ can be enhanced, when the level 
A in the 5 well coincides with one of the discrete levels 
in the well at i > xq. The Wigner resonant tunneling 
occurs in this case |12| . 

A scenario becomes different, if we return to the full 
potential energy l[2"]). After coming from under the bar- 
rier the electron participates in cyclotron motion and the 
Lorentz drift similar to a classical particle. The Lorentz 
drift in the potential (woj/ 4 /a 4 — £qx) brings the electron 
to the infinite x. This happens, according to electrody- 
namics, as a result of multiple reflections from the "wall" 
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uoy 4 /a, 4 [I3|- In that case Wigner tunneling does not 
occur since there are no discrete levels to the right of 
the barrier. This situation corresponds to a decay of the 
metastable state in the S well since after tunneling the 
electron goes to the infinite x. 

Whereas the potential y strongly disturbs the motion 
at the classical region, the underbarrier motion along the 
direction y = is hardly violated since reflections from 
that potential under the barrier are less important. It 
is shown in Sec. |V] For this reason, for an underbarrier 
motion in the full potential (f5]) one can use the approach, 
when the potential y 4 is dropped down. Accordingly, 
in the semiclassical approximation the tunneling rate is 
given by the WKB formula 



w ~ exp 




p 2 In 



(22) 



which follows from Eq. (|18p as \ip(xo,y)\ 2 . 

For zero magnetic field, h — 0, the tunneling rate (j22|) 
does not depend on A and coincides with the conventional 
WKB expression ©. For B = 15 and A = — 1 the tunnel- 
ing rate in Fig. [2] drops down with the magnetic field, as 
usually, since the metastable state in the quantum wire 
has zero momentum (|16|) perpendicular to tunneling. For 
B = 15 and A > — 1 such a momentum is finite resulting 
in cyclotron enhancement of tunneling. In Fig. [2] we plot 
the tunneling rate for A = —0.594. This value is chosen 
to have a link to numerical data discussed below. 



IV. CLASSICAL TRAJECTORIES 

For the potential analytical solutions of the 
Schrodinger equation ([7]) and the Hamilton- Jacobi equa- 
tion (TTUJ) do not exist. Nevertheless, to calculate a tun- 
neling rate with the exponential accuracy it is not neces- 
sary to solve the Hamilton- Jacobi equation (fit))) in the full 
{x, y} plane. It is sufficient to track a(x, y) along a clas- 
sical trajectory {x(r), j/(r)}, where r = —it is imaginary 
time since an underbarrier classical motion is impossible. 

The method of classical trajectories in imaginary time 
is well developed for static potential barrier with no mag- 
netic field [ill, Ell Hl| ■ In this case the main contri- 
bution to a tunneling rate comes from the extreme path 
linking two classically allowed regions. The real under- 
barrier path {x(r),y(r)} can be parametrized as a clas- 
sical trajectory in imaginary time. 

A situation with a magnetic field is substantially dif- 
ferent. For the potential ([2]) in a magnetic field the coor- 
dinate x(t) remains real in the underbarrier motion but 
y(r) = —vti{t) becomes imaginary. This means that the 
trajectory does not track the entire underbarrier path 
as for zero magnetic field. It just provides a "bypass" 
through the plane of complex y. So the method of clas- 
sical trajectories for tunneling in a magnetic field is non- 
trivial. 



One should note that the same situation takes place for 
tunneling across nonstationary barriers where imaginary 
time stands instead of imaginary y. Validity of the non- 
trivial trajectory method for the nonstationary barrier 
was numerically proved in Ref. 22]. The numerical re- 
sults of the present paper support the trajectory method 
also for a static magnetic field. 

We consider tunneling from the S well (r = tq) to 
the classically allowed region (r = 0). The ratio of the 
densities 



to 



l^(0),0)| ; 

hMo,o)| 2 



(23) 



can be identified with a tunneling rate. With the semi- 
classical accuracy, w ~ wt, where 

w T = exp (-2.BIm [tr(x(0), 0) - tr(0, 0)]) (24) 

is the tunneling exponent. The trajectory method allows 
to calculate the part 

A = 2Mm [a(x(0), 0) - cr(0, -ir)(iT ))] (25) 



of the total action (|24[) only. This part is expressed 
through the classical trajectory 



A Q = 2B 
drj 

-hx — 

OT 



(It 



4-A 

e 4 



See, for example, [l0|, The expression in the square 
brackets is the Lagrangian in terms of imaginary time. 
The coordinates x(t) and t)(t) in Eq. (|26|) are solutions 
of the classical equations of motion 

1 d 2 rj dx 



Id 2 



2 <9r 2 Or 



A 4 w 



with the conditions 
dx 



x(t ) = 0, 



77(0) - 0, 



dx 



0. (28) 



Since the trajectory terminates at the 5 well, this leads to 
the first condition ([28]) . The second condition (|28|) results 
from Eq. (TTTj) . The third condition (|28p corresponds to 
the physical exit point y = 0. The fourth condition l|28p 
is a property of an exit point where the electrons stops to 
move in the tunneling direction. Note, that the physical 
velocity perpendicular to the tunneling direction at the 
exit point from under the barrier is not zero. In the 
dimcnsionlcss units used it is 



dy_ 

dt 



(h) 



(29) 



A solution of Eqs. (|2T|) corresponds to the total energy 
A = -^l^| 2 + -^V-, + ^. (30) 



dx 



The relation (f3"01) is the fifth condition to the equations 
(|27p which can be satisfied by a proper choice of tq. 



5 




FIG. 3: The modified potential with the infinite wall at x — R. 
The discrete energy levels A„ are associated with the 5 well. 
The levels A^ are connected with the triangular well. 



o 



o> -8 




Total action 



The trajectory terminates at the unphysical (complex) 
point x — 0, y — —irj(To). One should connect this point 
with a physical one, for example, x — 0, y = 0. In other 
words, one should add to Aq the action 



A 1 = 2Blm [a(0, -ir)(i,T )) - cr(0, 0)] . 



(31) 



to complete the tunneling exponent (|24[) . One can find 
the action (|3"Tj) by a direct solution of the Hamilton- 
Jacobi equation (fT0| . Since the condition (fTT|) holds at 
all y, it follows from (fTTTj) that 



1 + A. 



The integration results in 



Ai = 2Be(l + A) 3 / 4 / 



»?(to) 


/ft? 




.£(l + A) 1 /4_ 




J 



(32) 



(33) 



FIG. 4: (Color online) Decimal logarithm of the electron den- 
sity at B = 15 calculated by the method of Sec. [V] The 
dashed curve: Zero magnetic field, h = 0. The solid curve: 
The finite magnetic field, h = 1.1, and the eigenvalue at the 
wire An = —0.594. The two dots mark the classical exit points 
calculated from Eqs. (fTSjl and (fT9|) . 



which follows from Eqs. ([2T]) . if to drop down the nonlin- 
ear part in 77. The total underbarrier time is determined 
by the relation 



1. 



p sinh (2/jto) 

One should also put f(z) = z in Eq. 
time To, as follows from Eqs. (|36f and 



1 + ^/T+p 1 



(37) 

. At the terminal 
371), 



(38) 



where 



NUMERICAL CALCULATIONS BY THE 

X - Y NET 



dv\f\ 



(34) 



The sign term in Eq. (|33[) accounts for a correct sign of 
the square root just to match the solution at the point 
x = 0, y — —ir){vo). So the tunneling rate is given by 



cxp(-A - Ai) , 



(35) 



where the parts of the total action are determined by 
Eqs. flU and ([33]). 

It is not difficult to check that, if to formally drop down 
the part y 4 in the potential (J2), the trajectory formalism 
(l35|) gives the same result (|22|) . For this purpose one 
has to write down the explicit expression for the classical 
trajectory 



sinh (2/ir) t 
h sinh (2h to) ft' 



(36) 



For numerical studies of the problem we consider the 
Schrodinger equation ([7]) at positive x with the boundary 
condition (|12p . Since in numerical calculations an infi- 
nite x is impossible we modify the potential U(x,y) @ 
adding the infinite potential wall at x = R. The wall is 
accounted for by the condition ip(R, y) — 0. In Fig. [3] the 
modified potential Um{x,Q) is shown, where the electron 
is completely localized resulting in discrete eigenvalues 
of energy A. In the S well, due to the potential y 4 /s 4 , 
there are discrete levels A„. With no y 4 potential the set 
A n becomes continuous. The levels associated with the 
triangular well are AW. 

For the potential U(x,y) $2$ an electron, initially lo- 
calized at the 5 well, can go to the infinite x providing 
a decay of this metastable state during the typical time 
to ~ 1/wt, where wt is the tunneling exponent (|24p. 

For a state, initially localized at the 5 well of the po- 
tential Uu{x,y), the same type of time to ~ 1/wt is 
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FIG. 5: (Color online) The three-dimensional plot of the deci- 
mal logarithm of the electron density for h — 1.1, A = —0.594, 
and B = 15. 



required to spread out to the triangular well when X n 
is not equal to A'^', that is we are away from Wigner 
resonances. The ratio of steady wave functions at the 
triangular well and at the S well is exponentially small. 
If to avoid Wigner resonances, it is determined by the 
same exponent as the tunneling decay rate for the case 
of an infinite R. This justifies the use of Um(x, y). 

When an energy level, occupied by the electron at the 5 
well, coincides with one at the triangular well the Wigner 
resonant tunneling occurs. In this situation the ratio ([23]) 
of steady wave functions at the triangular well and at the 
5 well is not exponentially small. Nevertheless, the time 
to fill out the initially empty triangular well still remains 
exponentially long, to ~ 1/ ^/wt- 

We numerically calculated ip(x,y) for the potential 
Um{x, y) using the discrete two-dimensional net in {x, y} 
space. A step was chosen from a condition of a bet- 
ter convergence. The method used works well for hyper- 
bolic system and does not result in unphysical exponents. 
We started with the asymptotic ^ ~ exp (— By 3 /3e) at 
y — > oo and chosen eigenvalues of A from the condition of 
d\ip\/dy = at y = 0. In the calculations we put B = 15, 
£ = 1, and R = 2. The calculations are performed at the 
interval < x < 2 with the boundary condition (fT2|) 
at x = 0. A distance between energy levels in Fig . [3] 
can be estimated from the semiclassical approach [12[ as 
(A„+i - A„) ~ (A( w+1 ) - AW) ~ 0.1 - 0.2. This is con- 
firmed by the numerical calculations. 

The results are figured out in Fig. [U 

At zero magnetic field, h = 0, the variables are sepa- 
rated and ip(x, 0) is a solution of the effective Schrodinger 
equation (Thi]) with the potential ((TSJ) , where one has to 
add the wall at x = 2 and formally put fe = h = and 
A ~ —1. The wave function is shown in Fig. U] by the 
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FIG. 6: (Color online) Decimal logarithm of the maximal 
value of the wave function after an exit from under the bar- 
rier versus magnetic field h. The maxima relate to Wigner 
resonances, A n = A' . A level in the wire from the interval 
0.55 < A n < 0.6 is adjusted to the variable h. The dashed 
curve is lg«j from Fig. [2] with A = —0.594 for comparison. 
It corresponds to the electron density after the barrier when 
R = oo (no Wigner resonances) . 



dashed curve. Due to reflections from the wall at x = 2 
the wave function oscillates at the interval 1 < x < 2. 
At the classical exit point x = 1 the wave function, cal- 
culated by Eq. (fl~8)) . is marked by the dot on the dashed 
curve in Fig. |4]where |t/>(1,0)| 2 ~2x 10~ 9 . 

At the finite magnetic field, h = 1.1, the numerically 
calculated wave function is plotted by the solid curve 
in Fig. |U The chosen eigenvalue A = —0.594 is one of 
the discrete levels A n in Fig. [3] The wave function is 
extended over the whole interval < x < 2 due to re- 
flections from the potential y 4 and the wall at x — 2 in 
the magnetic field. A classical exit point {ir(0),0} and 
4>(x(0),Q) can be calculated on the basis of the trajec- 
tory method of Sec. IIVI The results hardly differ from 
the case of Sec. UTTj when the potential y 4 was dropped 
down. According to Eq. (|3"5)) , 77(70) — 0.263 which leads 
to the very small potential y 4 ~ 0.005. It is obvious and 
also can be shown that an influence of this potential on 
the classical trajectory is small. We emphasize that in 
our case the potential y 4 hardly influences the underbar- 
rier motion only. After an exit from under the barrier 
the electron is strongly reflected by that potential. 

Neglecting the part y 4 , one can easily perform trajec- 
tory calculations. This results in the exit point (I19p . 
x = x(0) ~ 0.594, and the electron density \ip(xo, 0)| 2 ~ 
2.0 x 10~ 5 marked as the dot on the solid curve in Fig. [4] 
We see that the analytical and numerical results are close 
to each other for both zero and the finite magnetic fields. 

The three-dimensional plot of \ip(x, y)\ 2 is shown in 
Fig. [5] Due to reflections from the potential y 4 and the 
wall at x — 2, the motion covers the whole region apart 
the barrier. We also performed numerical calculations 
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for the potential y 16 instead of y just to be closer to a 
situation of an infinite restricting wall. Qualitatively, the 
wave function looks the same filling out the whole {x, y} 
space to the right of the barrier. 

To get tunneling to the triangular well more adequate 
to tunneling to an infinite region, we should avoid Wigncr 
resonances, as pointed out above. In Fig. [5] the maximal 
value of \ip(x, 0)|^ laa , after the exit point is plotted versus 
h. As one can see, the structure of Wigner resonances is 
very pronounced. A shape of the wave functions in the 
triangular well close to Wigner peaks was checked. We 
found no deviation from eigenfunctions at the triangular 
well, when the S well at x — was absent. This proves 
the Wigner nature of the peaks in Fig. [6] We choose 
h = 1.1 to get away of the Wigner resonances as clear 
from Fig. [6] 



VI. NUMERICAL CALCULATIONS BY THE 
MATRIX FORMALISM 



™ -4 



V ' 

VS. a 


1 ■ 


■ 

- 


X \ 
















i J 



0.0 0.5 1.0 1.5 2.0 

x 

FIG. 7: (Color online) Decimal logarithm of the electron 
density at B = 15 calculated by the method of Sec. IVII 
The dashed curve: Zero magnetic field, h = 0. The solid 
curve: The finite magnetic field, h — 1.1, and the eigen- 
value of the wire state A n = —0.586. The two dots mark 
the classical exit points |i/>(l,0)| 2 ~ 2 x 10~ 9 for h = 
and 1-0(0.586, 0)| 2 ~ 2.3 x 10" 5 for h = 1.1 calculated from 
Eqs. (fig )) and l[T9 ]) , 



In Sec. |V] we performed the numerical calculations at 
the interval < x < 2 using the boundary condition (TT2"]) 
at x — 0. In this section we demonstrate the numerical 
formalism, which does not explore the boundary condi- 
tion (fT2|) . The calculations are performed at the interval 
—2 < x < 2, where the potential Um(%,u) is symmetri- 
cally continued from the region < x for negative x. The 
6 function is approximated by a narrow well of the width 
of one computational step. The same method to model a 
5 well was used in Ref. [22| . For positive x in the semiclas- 
sical approach (a large B) the results for the symmetric 
potential chosen and for the potential Um(x, y) should be 
close to each other. This is confirmed by the numerical 
calculations. 

In this section we use a matrix formalism. We start up 
with the Schrodinger equation 



1 9 h 
thy 

B ox 



ip - 



— - M - 15 5 ( x ) 
e B 



1 d 2 ip 
B^'dy 2 



(39) 



equations 



1 8 2 F n 2ih ^ dF m 



B 2 dx 2 5 3 / 2 



Ox 



m=0 x 

1 + 2n 
B 



\x\--8(x)-X 



F n = 0, (42) 



where the matrices a™, /3™, and 7™ are symmetric and 

dzip m (z)z 4 ip n (z). (43) 



The matrix /3™ (7™) is determined by the relation analo- 
gous to (|4"3"|) but with the substitution z 4 — > z (z 4 — » z 2 ). 
The non-zero elements are 

_ 4 _ y/n(n-l)(n-2)(n-3) 



with the conditions ip(±2,y) — 0. Eq. ([3^)1 differs from 
the form |(7J) by the gauge. The wave function can be 
written as the expansion 



„ _ 3(2n 2 + 2n + 1) n _ 2 _ (2n - 1) ^n(n - 1) 



i>(%,V) = ^2F n (x)(p n (yVB), 



(40) 



n=0 



where ip n (z) is a normalized eigenfunction of a harmonic 
oscillator, 



d 2 ip n 
dz 2 



+ z 2 ip n = (1 + 2n)ip T , 



(41) 



One can easily show that the functions F n (x) satisfy the 



if 



yjn(n - 1) 



l + 2n 



(44) 



and also a™ +4 , a™ +2 , /3^ +1 , and 7^ +2 obtained as sym- 
metric combinations from (|4"4"|) with the shifts of n. 

The above matrix formalism is used to numerically cal- 
culate the wave function (|4"D|) by a solution of the eigen- 
value problem (|4"2")1 . The calculations were performed by 
means of the discretization of the interval — 2 < x < 2 
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using 513 points. To check the scheme we made addi- 
tional runs for the interval — 3 < x < 1 with the equal 
number of discrete points. The both methods resulted in 
the same eigenvalues and eigenf unctions. In all the cal- 
culations we considered nine terms in the expansion (|40|) 
of the wave function. 

The results are displayed in Fig. [7| showing the states 
generic with ones in Fig. 21 The difference in A n for 
the both figures is 0.008 which is substantially less than 
(A n +i — X n ) ~ 0.2. 

Comparing Figs. [7] and one can conclude that two 
different numerical methods lead to very similar results. 
The numerical data of the both methods are also in a 
good agreement with the trajectory results shown by the 
dots in Figs. Hand El 



VII. DISCUSSIONS 

A common point of view is that due to intrinsic un- 
derbarrier mechanisms a magnetic field always reduces 
a probability of tunneling decay of a metastable state. 
In other words, it never chances to get tunneling rate en- 
hanced by a magnetic field for many orders of magnitude. 
This point of view is supported by a general arguments 
that in the classical electrodynamics a motion ahead is 
prevented by cyclotron rotation and the same happens 
in quantum mechanics. The famous example is Landau 
states. The increased localization of a wave function in a 
magnetic field can also be explained by an additional po- 
tential, proportional to squared velocity, created by the 
magnetic field. Any hand waving explanation leads to a 
conclusion of tunneling reduction. We propose a mecha- 
nism which contrasts to that point of view. 

Cyclotron enhancement of tunneling is clear for under- 
standing. It occurs due to a rotation under the barrier 
of a velocity vector of the electron. Under the cyclotron 
rotation it becomes more aligned with the tunneling di- 
rection, which increases the tunneling rate. A necessary 
condition for the phenomenon is a finite velocity (JTTJ) 
in the decayed state of the quantum wire. The velocity 
should be directed along the vector £o x H , which is par- 
allel to the wire. In the case of the opposite direction the 
magnetic field reduces the tunneling rate. For cyclotron 
enhancement of tunneling it is important that a sample, 
from where tunneling occurs, should be extended in the 
direction perpendicular to tunneling. The effect is absent 
in tunneling from a small quantum dot. 

The scenario of decay of the metastable state in the 
magnetic field consists of two parts. First, the electron 
moves under the barrier rotating its velocity. On this 
step, a role of reflections from the restricting walls is mi- 
nor. Second, after exit from under the barrier the elec- 
tron goes to infinity performing multiple reflections from 
the restricting walls. 

The semiclassical approach used is confirmed by two 
different numerical methods. A domain of parameters 
for the phenomenon can be roughly estimated from the 



condition that the Lorentz force, caused by the velocity 
in the wire, is of the order of the potential force £q. The 
parameter h is 



h = B- 



2uo 



, iJ(Tesla) 
0.58 x IQ- 4 B- v 



M (eV) 



(45) 



Since we are interested in h ~ 1 and the semiclassical 
parameter B is large, the Landau splitting ftw c is always 
less then uq which defines a scale of the energy levels in 
the quantum wire. 

To be specific, let us consider different regimes of tun- 
neling. 

Strong effect on tunneling (h = 2). We take B = 40 
and A = —0.3. Then at zero magnetic field the tunneling 
rate (HH) is 6.9 x 10~ 24 . At h = 2.0 the tunneling rate Igg) 
becomes 1.1 x 10~ 7 . The enhancement occurs under the 
condition uo(eV) ~ 1.16 x 10 _3 -ff (Tesla). For reasonable 
values of the magnetic field, H ~ 10 Tesla, the quantum 
wire or the thing film should be "soft" in the sense of a 
not large u ~ 0.01 eV. See Refs. QJ, EE HI- The exit 
point from under the barrier is estimated as xq ~ 250A 
and the electric field is £q ~ 1.3 X 10 3 eV/cm. 

Weak effect on tunneling (h — 0.2). We take B = 15.8 
and A = —0.3. Then at zero magnetic field the tunneling 
rate (HH) is 7.1 x 10~ 10 . At h = 0.2 the tunneling rate (|2"2"|) 
becomes 1.0 x 10~ 7 . The enhancement occurs under the 
condition uo(eV) ~ 1.16 x 10~ 2 i? (Tesla). For reasonable 
values of the magnetic field, H ~ 10 Tesla, one can choose 
ito ~ 0.1 eV. The exit point from under the barrier is 
estimated as xq ~ 75A and the electric field is £o ~ 
1.0 x 10 5 eV/cm. 

One can conclude that an experimental observation of 
cyclotron enhancement of tunneling is possible. An ex- 
perimental arrangement can be a quantum wire embed- 
ded into a two-dimensional electron system. The electric 
field is parallel to the two dimensional system and is per- 
pendicular to the wire. The magnetic field is perpendic- 
ular to the wire and to the electric field. One can use a 
dielectric wall parallel to the tunneling direction just to 
get electrons moved away after tunneling along the wall 
due to multiple reflections from it in the magnetic field. 

Another way of an experimental observation is to use a 
thin film, with a quantized motion inside, perpendicular 
to the electric field. 



VIII. CONCLUSIONS 

A magnetic field can enhance a tunneling decay rate of 
a metastable state for many orders of magnitude. This 
happens due to intrinsic underbarrier mechanisms of a 
cyclotron rotation of an electron velocity. It becomes 
more aligned with the tunneling direction resulting in 
the enhancement of the tunneling rate. The effect can be 
observed experimentally for tunneling from a quantum 
wire or a thing film across a barrier created by an applied 
electric field. 
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